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I i Abstract 

o 

5— ( ■ We present near-infrared spectroscopic observations of star-forming galaxies at z 1.4 with EMOS on 

C/3 \ the Subaru Telescope. We observed K-band selected galaxies in the SXDS/UDS fields with K < 23.9 mag, 

1.2 < Zph < 1.6, > IO^-^Mq, and expected F{Ha) > lO"^*^ erg s^^ cm"^ objects in the sample have 
significant detections of Ha. Eor these objects, excluding possible AGNs identified from the BPT diagram, 
gas-phase metallicities are obtained from [Nil] /Ha line ratio. The sample is split into three stellar mass 
' bins, and the spectra are stacked in each stellar mass bin. The mass-metallicity relation obtained at z ~ 1.4 

■"nI" , is located between those at z ^ 0.8 and z ~ 2.2. We constrain an intrinsic scatter to be '-^0.1 dcx or larger in 

' the mass-metallicity relation at z 1.4; the scatter may be larger at higher redshifts. We found trends that 

the deviation from the mass-metallicity relation depends on the SER and the half light radius: Galaxies 
with higher SER and larger half light radii show lower metallicities at a given stellar mass. One possible 
scenario for the trends is the infall of pristine gas accreted from IGM or through merger events. Our data 
points show larger scatter than the fundamental metallicity relation (EMR) at z ~ 0.1 and the average 
metallicities slightly deviate from the EMR. The compilation of the mass-metallicity relations at z ~ 3 to 
z 0.1 shows that they evolve smoothly from z ~ 3 to z '-^ without changing the shape so much except 
for the massive part at z 0. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: abundances 



1. Introduction lished the mass-metallicity relation at z ~ 0.1 with a large 

sample (- 53, 000) of Sloan Digital Sky Survey (SDSS) 

The gas-phase metallicity (hereafter, metallicity) is a galaxies. Tracing the cosmological evolution of the mass- 
fundamental parameter for the understanding of forma- metallicity relation is indispensable to reveal how galaxies 
tion and evolution of galaxies, because it traces past star have been evolving. At z ^ 1, the mass-metallicity rela- 
formation activity; metals in gas are produced in stars tions and the downsizing- like evolution to z ~ 0.1 were 
formed and returned into the inter-stellar medium (ISM) found; the less evolution can be seen in the massive part 
of galaxies. The metallicity is also affected by inflow and (e.g., Savaglio et al. 2005; Zahid, Kewley, & Bresolin 
outflow of gas. Thus, by investigating the metallicity and 2011). However, the anti-downsizing-like evolution of the 
its cosmological evolution, the star formation history of mass-metallicity relations were also presented at similar 
galaxies together with the gas inflow and outflow can be redshifts (Lamareille et al. 2009; Pcrez-Montero et al. 
constrained. 2009). The mass-metallicity relations were obtained at 

The presence of a correlation between stellar mass and z ^ 2 by using ~ 90 galaxies (Erb et al. 2006) and at z '--^ 3 

metallicity (hereafter, mass-metallicity relation) in the lo- by ~ 20 galaxies (Maiolino et al. 2008; Mannucci et al. 

cal universe is well known. Tremonti et al. (2004) estab- 2009) and the smooth evolution of the mass-metallicity 
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relation from z ~ 3 to 2 ^ is suggested. However, at 
z ^ 2, higher metallicities at a given stellar mass have 
been reported (Hayashi et al. 2009; Yoshikawa et al. 
2010; Onodera et al. 2010) than those found by Erb et al. 
(2006). The discrepancy at z ^ 2 may be partly due to 
the small size of their samples (10 — 20). Since the red- 
shift of z = 1 — 2 is close to the peak epoch in the cosmic 
star-formation history, establishing the mass-metallicity 
relation at this redshift is very important. 

Although the metallicity correlates well with stellar 
mass, there is a scatter in the relation. Its origin may pro- 
vide a clue to understand the process of the chemical en- 
richment and may be related to how the mass-metallicity 
relation evolves. At z '-^ 0.1, Tremonti et al. (2004) found 
that the mass-metallicity relation has an intrinsic scat- 
ter of ^ 0.1 and the scatter is larger at the lower stellar 
mass. At a given stellar mass, galaxies with higher surface 
stellar mass density tends to show higher metallicity, sug- 
gesting that they transformed more gas into stars raising 
the metallicity. Ellison et al. (2008a) showed that galax- 
ies with larger specific star-formation rate (SFR) and size 
show lower metallicity at a given stellar mass at 2 ^ 0.1. 
More recently, Mannucci et al. (2010) and Yates et al. 
(2011) found that galaxies with larger SFRs tend to show 
lower metallicities at z 0.1 and the scatter around the 
mass-metallicity relation is reduced very much (i.e., the 
fundamental metallicity relation) by introducing SFR as 
the second parameter affecting the metallicity. The scat- 
ter of the mass-metallicity relation, however, is still un- 
clear at higher redshifts because of the limited sample size. 
In order to investigate these trends further, including the 
dependency of other parameters, at higher redshift, a large 
near-infrared spectroscopic survey is necessary, which is 
very time consuming and hard to achieve. 

The Fiber Multi Object Spectrograph (FMOS) on the 
Subaru Telescope enables us to survey a large spectro- 
scopic sample at high redshifts. FMOS is a near infrared 
(NIR) fiber multi-spectrograph (Kimura et al. 2010). The 
fiber positioner "Echidna" at the prime focus feeds 400 
fibers in a FoV of 30' diameter. Two NIR spectrographs 
(IRSl and IRS2) with an OH airglow suppression system 
are capable of obtaining both low resolution (i? ~ 650) 
and high resolution {R ~ 3000) spectra in the wavelength 
range of 0.9 — 1.8/im. Taking advantage of the multiplicity 
available with FMOS, we are conducting a large spectro- 
scopic survey for star- forming galaxies at z = 1 — 2. This 
redshift range is of much interest, because it is near or 
shortly after the peak epoch in the cosmic star formation 
history (e.g., Hopkins & Beacom 2006) and thus impor- 
tant for the history of the chemical evolution of galax- 
ies. It is also suitable for FMOS, because in the redshift 
range [Nll]AA6548,6584, Ha, [O III] AA4959,5007, and H^ 
lines are located in the wavelength range observable with 
FMOS. In this paper we present the first results from 
this survey as to metallicity measurements of galaxies at 
z ~ 1.4. 

Throughout this paper, we adopt the concordance cos- 
mology, (r^A/ , r^A , /i) = (0.3, 0.7, 0.7). All magnitudes 
are in the AB system (Oke & Gunn 1983). 



2. Sample and Observations 

2.1. Photometric Data in the SXDS/UDS Fields 

The sample is selected from a galaxy catalog in the 
SXDS/UDS fields, where deep multi- wavelength data are 
available. The optical ^MfearM/Suprime-Cam images (B, 
V, Rc, and z') are available from the Subaru-XMM 
Deep Survey (SXDS; Furusawa et al. 2008) with the orig- 
inal pixel scale of 0"202 pix~^. The limiting magnitudes 
are 27.7, 27.2, 27.1, 27.0, and 26.0 mag for B, V, Rc, 
and z'-band, respectively (5cr, 2"0 diameter). The NIR 
UKIRT /WFCAM images (J, H, and K^) are taken from 
DR8 version of the UKIDSS Uhra Deep Survey (UDS; 
Lawrence et al. 2007) with the original pixel scale of 0."268 
pix~^. The limiting magnitudes are 24.9, 24.2, and 24.6 
mag for J, H, and Ks, respectively (5cr, 2"0 diameter). 
The Spitzer/IRAC images (3.6/Ltm, 4.5/xm, 5.8/xm, and 
8.0/Ltm) are taken from the Spitzer public legacy survey of 
the UKIDSS Uhra Deep Survey (SpUDS; PI: J. Dunlop, in 
preparation) with the original pixel scale of 0."600 pix~^. 
The limiting magnitudes are 24.8, 24.5, 22.8, and 22.7 mag 
for chl (3.6/^m), ch2 {4.5^im), ch3 (5.8Aim), ch4 (8.0/im), 
respectively (5cr, 2."4 diameter). We use an overlapped 
area of these data. The effective area is ^ 0.67 deg^. 

The NIR and MIR images are aligned to the optical 
images with pixel scale of 0."202 pix"-'^. Object detection 
and photometry in each image is carried out with respect 
to the reference i^T^-band image by using the double im- 
age mode of SExtractor (Bertin &; Arnouts 1996). For 
the optical and NIR images, aperture photometry with 
an aperture of 2"0 diameter is carried out after apply- 
ing the PSF matching. The total magnitude is calculated 
from the aperture magnitude scaled to the SExtractor's 
MAG_AUTO in X^-band. For the IRAC images, the to- 
tal magnitude is calculated from the aperture magnitude 
with 2."4 diameter applying aperture corrections. The cor- 
rection factors are determined by generating artificial ob- 
jects with various intrinsic sizes on the IRAC images and 
recovering the aperture magnitudes. The uncertainties of 
the aperture correction in the above simulations are 0.05, 
0.05, 0.09, and 0.08 mag for chl, ch2, ch3, and ch4, re- 
spectively. 

2.2. Photometric Redshift Sample 

For each source, the photometric redshift (photo-z) 
is determined by using Hyperz (Bolzonella et al. 2000), 
where model templates are fitted to the observed Spectral 
Energy Distributions (SEDs) from optical to MIR by a 
standard minimization. We use a standard set with 
synthetic spectra with 51 age grids from 0.3 Myr to 20 
Gyr and 21 reddening {Ay) grids from 0.0 to 4.0 with 
AAv = 0.2 mag, and the CWW (Coleman et al. 1980) 
templates as well. The range of the photo-z is set to be 
0.0 < Zph < 6.0 with a step of Az=0.03. The obtained 
photo-z is compared with the available spectroscopic red- 
shift (spec-z) in the SXDS/UDS (Simpson et al 2011, in 
preparation; Akiyama et al. 2011, in preparation; Small 
et al. 2008). Although some catastrophic outliers exist, 
the photo-z generally agrees with the spec-z with uncer- 
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Fig. 1. Our targets on the Af»— SFR diagram. Left: Our primary sample with Ks < 23.9 mag, 1.2 < Zp^ < 1.6, and 
Alt > 10^'^Mq is indicated by gray dots. Among the primary sample, the secondary sample with expected F(Hq!)>1x 10^^" 
erg s~^ cm~^ is indicated by blue open circles. Right: Objects observed with FMOS are plotted as circles. While 
objects whose Ho lines are detected with SN of >3 are indicated by red filled circles, objects whose Hq is not de- 
tected with SN of >3 are indicated by blue open circles. The secondary sample are also plotted by gray dots. 



tainty of aj{l + z) = 0.02 - 0.03 at z = 1 - 2. 

In order to derive stellar masses of these galaxies, we 
perform SED fitting by using the SEDfit (Sawicki 2011). 
For model SEDs, we use the BC03 population synthesis 
code (Bruzual & Chariot 2003) with Padova 1994 evolu- 
tionary tracks. The Salpeter IMF (Salpeter 1955) with 
the mass range of 0.1 — 100 M© and the mctallicity of 1.0 
Zq are assumed. The extinction curve of Calzctti et al. 
(2000) is also assumed. We fix redshifts to those deter- 
mined by the Hyperz in order to reduce the number of 
free parameters. We generate model SEDs with 25 age 
grids from 0.2 Myr to 15.8 Gyr, 5 star formation histories 
(constant SF, r = IGyr, r = lOOMyr, r = lOMyr, and in- 
stantaneous burst) grids, and 21 E{B — V) grids from 0.0 
to 1.0 with a step of AE(B ~V)== 0.05 mag. Wc then fit 
the model SEDs to the observed SEDs with the standard 
minimization. 

The color excesses are estimated from the rest-frame 
UV colors. It is known that the rest-frame UV color 
strongly correlates with E{B — V) (Mcurcr et al. 1999). 
For objects at z 2, Daddi et al. (2004); Daddi et al. 
(2007) presented the relation between B — z' color and 
E{B — V) by calibrating with the results from the SED 
fitting. For the redshift range of our sample, B — z' color is 
not a good indicator of the extinction because the z'-band 
covers beyond the Balmer/4000Abreak. For 1.2 < z< 1.6, 
we calibrate the relation between the B ~ i' color and 
E{B - V) by using the rcsuhs of the SED fitting. The 
resulting relation is E{B -V) = 0.31(B - i') + 0.02. AU 
the color excess of our sample is derived from the B — i' 
color by using this relation. As we mention in Section 3, 
both Ha and H/3 lines are detected in some galaxies. For 
these objects, the color excess can be derived based on the 



Balmer decrement, by assuming the intrinsic Ha/Hp ratio 
of 2.86 and the Calzetti extinction law. The average re- 
sult is consistent with the color excess from the rest-frame 
UV colors within the error, if the extinction for the neb- 
ular emission is assumed to be larger than that for stellar 
continuum as we describe below. 

The SFR is derived from the rest-frame UV luminos- 
ity density by using the conversion by Kennicutt (1998). 
The rest-frame UV luminosity density is calculated from 
the i3-band magnitude. The intrinsic SFR is derived by 
correcting for the extinction with the E{B — V) derived 
above. The Calzetti extinction curve is also assumed. For 
some objects, we cross-check with the SFRs from MIPS 
fluxes by using archived images of the SpUDS. We use 
the conversion from the rest-frame Ls/jm to the SFR by 
Daddi et al. (2007). The resuhing SFR from IR lumi- 
nosity roughly agrees with that from the rest-frame UV 
luminosity density. 

The expected Ha flux is calculated from the intrinsic 
SFR and the E{B — V) described above. Since it is sug- 
gested that the extinction is significantly larger for the 
ionized gas than for the stellar component (Calzetti et al. 
2000), we convert the obtained E{B — V) to that for the 
ionized gas by using a prescription by Cid Fernandes et al. 
(2005). For the conversion from the SFR to the Ha lumi- 
nosity, we use the relation by Kennicutt (1998). 

The sample is cross-correlated with X-ray sources 
(Ueda et al. 2008) and objects which are cross-matched 
within the error circle of the X-ray source arc excluded 
from the sample. Accordingly, the X-ray bright AGNs 
(-'^x(2-iokcV)^"'^^^'^ ^''S ^""^^ excluded from our sam- 
ple. 
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2.3. Sample 

For the spectroscopic observations with FMOS, the pri- 
mary sample is constructed from the K-selected cata- 
logue described above with the following selection criteria: 
Ks < 23.9 mag, 1.2 < Zph < 1.6, and > IO^-^Mq. This 
primary sample size is very large, and many of the sample 
galaxies may show very faint Ha emission lines. Thus in 
order to make an efhcient survey, and as a first step, we 
selected objects for which expected Ha flux is larger than 
1.0x10"^^ erg s~^ cm~^ as the secondary sample. This 
selection could introduce some bias for the sample, hence 
the effects will be considered later based on the results 
obtained. 

In Figure 1, intrinsic SFRs derived from the rest- 
frame UV are plotted against the stellar mass (hereafter 
M^, — SFR diagram) for the primary sample. The SFR 
correlates with the stellar mass. On the Af*— SFR dia- 
gram, the secondary sample is also shown with blue open 
circles. The distribution is slightly shifted toward larger 
SFR at a given stellar mass particularly in less massive 
part. This is due to the fact that we selected objects with 
larger expected Ha flux and thus larger intrinsic SFR. We 
will discuss the effect of the selection bias on the metallic- 
ity in Section 3.9. Actually observed objects are randomly 
selected from the secondary sample. The right panel of 
Figure 1 shows the observed objects among the secondary 
sample. 

Some target objects are relatively massive (M* ^ 5 x 
10^° Mq) and have high intrinsic SFR (> lOOOMeyr"!) 
as seen in Figure 1. We cross-correlated to the MIPS 
sources, and we found that the detection rate in MIPS 
is significantly lower for these objects than for other ob- 
jects. The estimation of these high SFRs may be incorrect 
partly because the E{B — V) is overestimated due to the 
dust/age degeneracy. We included these objects in our 
targeted sample and found that the actual detection rate 
of Ha emission for the galaxies in this A/* — SFR region 
is significantly lower than that in the other region as seen 
in right panel of Figure 1. 

2.4. Observations and Data Reduction 

The observations were carried out with FMOS (Kimura 
et al. 2010) on the Subaru Telescope on September 27 — 28, 
November 11-13, 21-24, 28, and December 15-16, 2010 
(mostly the guaranteed time observing runs and partly en- 
gineering observing run for science verification, and open 
use observations). The observations were made with the 
Cross Beam Switch (CBS) mode in low resolution (LR) 
mode employing both IRSl and IRS2. The spectral reso- 
lution in the LR mode is i? ~ 500 at A ~ 1.10/im, R ~ 650 
at A ~ 1.30/xm, and R ~ 800 at A ~ 1.55/xto, which are 
measured from the Th-Ar lamp frames. The pixel scale in 
the LR-mode is 5 A. Our targets are allocated to fibers 
together with other scientific targets. Although weather 
conditions varied during the observing runs, the typical 
exposure time was ^ 3 hours for one FMOS field of view. 
The typical seeing size measured with the FMOS sky- 
camera during the observations was 0."9 in i?-band. 



The obtained data are reduced with the FMOS pipeline 
FIBRE-pac and detailed descriptions for the data reduc- 
tion arc presented by Iwamuro et al. (2011). Here we de- 
scribe an outline of the process. For a set of exposures, the 
A — B sky subtraction is carried out. Then the obtained 
2D spectra are corrected for the distortion and residual 
sky subtraction are also carried out. The obtained spectra 
are combined for the total exposures. Since the spectrum 
of one target is obtained by two fibers in the CBS mode, 
the final spectrum of one target is obtained by merging 
the CBS pair spectra. The wavelength calibration is done 
by using Th-Ar lamp frames. The uncertainty associated 
with the wavelength calibration is AA ~ 5A. The rela- 
tive flux calibration was done by using several F, G or 
K-type stars selected based on J — H and H — Kg colors 
and observed simultaneously with other scientific targets. 
The uncertainty of the relative flux calibration is ~ 10% 
at A = 1.0 — l.Sfim. The absolute flux is determined from 
the observed count rate by assuming the total throughput 
of the instrument, which is calibrated by using the moder- 
ately bright stars (~16 mag in J- or i?-band) in the best 
weather conditions ('^ O'.'l) in previous engineering obser- 
vations. The uncertainty of the absolute fiux calibration 
itself is estimated to be ^ 10%. However, there may be ad- 
ditional uncertainties (20 — 30%) that come from variation 
of the atmospheric transmission and the seeing conditions. 

In total, 317 objects were observed and for 71 objects 
the Ha emission line was detected with a signal-to-noise 
ratio (S/N) larger than 3, where the noise level is mea- 
sured from the continuum in a wavelength window of 
±0.1/im from the emission line. Detection rates are 26%, 
44%, 21%, and 9% at < 10^°Mq, IO^^Mq < < 
IO^O-^Mq, IO^O-^Mq < < 10"Mq, and > IO^Mq, 
respectively. Examples of the obtained spectra are pre- 
sented in Figure 2. Spectra in which the S/N of [N II] A6584 
is larger than 3.0, 1.5 < S/N < 3.0, S/N < 1.5 are shown 
in the left, middle, and right panels of the figure, respec- 
tively. 

3. Results 

3.1. Model Fitting to Emission Lines 

In order to measure a spec-z, a fiux of an emission line, 
and a line width, we made a model fit to the reduced one- 
dimensional emission line profile. This process is straight- 
forward in the usual spectroscopy. In our case, however, 
we need rather complicated processes. FMOS uses an OH 
suppression system to remove the strong OH airglow sky 
emission lines. The dispersed light is focused onto the 
mask mirror with a resolving power of i? ~ 3000 (~ 100 
km s~^) and is rejected by the OH masks. Then the spec- 
trum is anti-dispersed by VPH grating to i? 650 (~ 450 
km s~^) and focused onto the detector (see Kimura et al. 
(2010) for the optical design of FMOS). Since the reduced 
spectra are divided by a spectrum of the calibration star 
in the relative fiux calibration process of the pipeline, the 
effect of the OH-mask on the continuum light is corrected 
in the final reduced spectra. The correction, however, is 
not appropriate for emission lines if the line width is rel- 
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Fig. 2. Examples of the observed spectra {red solid lines) with the best-fitted models {blue dotted lines). In the 
left 5 panels, objects whose S/Ns of [Nll]A6584 lines are larger than 3 are plotted. In the middle panels, ob- 
jects with 1.5 < S/N < 3.0, and in the right panels, objects with S/N < 1.5 are plotted. Vertical dot- 
ted lines indicate the positions of [Nll]A6548, Ha, and [Nll]A6584, from left to right. The positions and the 
width of the OH-mask are indicated by hatched stripes. Continua are subtracted by third-order polynomial fitting. 
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atively small. Thus we present two different ways of the 
spectral fitting. One is fitting with model spectra taking 
into account the mask effects, and the other is fitting the 
observed spectra without taking the mask effects into con- 
sideration. In both cases, we subtracted the continuum 
by fitting with third-order polynomials in the wavelength 
range of 1.45— 1.70/xm. 

In the former method (methodl), we fit model spec- 
tra including the OH-mask effect to the observed spec- 
tra. The model spectra are generated as follows: An in- 
trinsic spectrum is modeled by a multi- Gaussian function 
with Ha and [Nll]AA6548,6584 hues. We assume that ra- 
tio of [Nll]A6584 and [Nll]A6548 is 3.0, the line widths 
are the same as that of Ha, and there is no systematic 
shift between Ha line and [Nil] lines. Each spectrum is 
parametrized by four parameters: redshift, amplitude of 
Ha, amplitude of [Nll]A6584, and line width. The spec- 
trum represented by these parameters is convolved to a 
resolution of ii ~ 3400, and then, the flux density in the 
regions (AA 8A) of the OH-masks is set to be zero. The 
obtained spectrum is convolved to a low-resolution spec- 
trum with i? - 800. 

In order to correct the continuum level affected by the 
OH-suppression, a flat spectrum in F^, without emission 
lines mimicking a star used in the flux calibration is pro- 
cessed in the same manner as described above. The correc- 
tion factors as a function of wavelength are derived from 
the resulting spectrum divided by the input flat spectrum. 
The model spectrum with emission lines is corrected by 
using this correction factor in order to compare with the 
observed spectrum. A schematic view of line profiles af- 
fected by the OH-mask is presented in Figure 3. 

The fitting process is carried out using mpfit pack- 
age (Markwardt 2009) which implements non-linear least 
square fitting with standard Levenberg-Marquardt algo- 
rithm. In the fitting procedure, we search for the mini- 
mum in the parameter space of the four parameters 
described above. The resulting maps of fits show that 
there are no pronounced discontinuities or strong degen- 
eracies between the parameters. 

For the latter method (methodE), we simply fit the ob- 
served spectrum by the multi-Gaussian profile. The fitting 
procedure is generally the same as methodl. The fitting 
models are multi-Gaussian functions, but without taking 
the mask loss corrections into consideration. 

In order to quantify the effects of the OH-masks and 
to evaluate the accuracy of the two methods, we generate 
artificial spectra with one emission line and no contin- 
uum, varying the line widths (50 - 1000 km s'^) and 
fluxes (1.0x10"^'' — 1.0 X 10"-^^ erg s~^ cm~^), which 
cover the ranges of widths and fluxes of most of the emis- 
sion lines detected with S/N > 1.5. Then we place an 
OH-mask virtually changing its wavelength randomly to 
simulate the amount of flux loss for an emission line at dif- 
ferent redshifts. The final LR spectrum is then made by 
degrading the spectral resolution down to 650. The 
random noise, which is equivalent to the average noise 
level actually observed in a typical exposure time, is added 
into each spectrum. The resulting artificial spectra are 
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Fig. 3. Schematic view of the OH-mask effects. The intrin- 
sic spectrum with At) of 200 km s~^ in the HR mode and 
that absorbed by the OH-mask arc indicated by a dotted line 
and a thin solid line, respectively. The absorbed spectrum 
degraded to the LR is indicated by a thin dashed line. The 
spectrum corrected for the mask effects by using continuum 
Hght is indicated by a thick solid line; this is the obtained 
spectrum after data reduction. The intrinsic spectrum de- 
graded to the LR without considering the mask effect is shown 
by a thick dashed line as a reference. See also text for details. 



fitted by two different methods [methodl and method2), 
and the obtained redshifts (wavelengths), line fluxes, and 
line widths are compared to the input values. We repeat 
this process 1000 times. The redshifts are recovered with 
an accuracy of ^ 5% in both methods. The line widths 
are underestimated up to ^ 50% in method2. While, in 
methodl^ the line widths can be recovered properly with 
an accuracy of < 20%. 

The recovery rate of the flux against the amount of the 
loss by the mask is presented in Figure 4. Our simulations 
show that method2 cannot properly recover the flux which 
is affected by the OH-mask, especially when the line width 
is small. For a line width of '^50 km s~^, up to ^90% of 
the line flux is damped by the OH-mask. In this case, the 
flux recovered by method2 is only ^20% of the intrinsic 
flux, and ^ 80% of the light is lost. The recovery of the flux 
improves if the line width is larger. For instance, the flux 
can be recovered within ^20% accuracy if the line width 
is ^600 km s~^. It is also interesting to note that the 
flux is overestimated by method2 in some cases, because 
the line profile is overcorrected by the correction for the 
continuum: If a large part a pixel covers an OH-mask 
and the peak of a narrow emission line just outside of the 
mask, the flux loss rate in the pixel for the continuum is 
larger than that for the emission line. 

Figure 4 shows that fluxes recovered by methodl agree 
well with the intrinsic fluxes for any line width. The uncer- 
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Fig. 4. The average ratio of recovered flux and the input flux against the amount of the loss by the OH-mask for intrinsic line 
widths of 50 km s~^ (top left)^ 100 km s""^ [top right)^ 200 km s~^ {middle left), 400 km s^^ (middle right), 600 km s""^ [bottom 
left), 1000 km s~^ (bottom right). Red lines and circles are the results by methodl, while the blue lines and triangles are those 
by method2. Error bars indicate Icr standard deviations of the simulated distributions. See text for details of the simulations. 
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taiiity of the correction is ^ 20% in most cases. If the loss 
is ^ 90%, however, the uncertainty is ^ 30%. Of course, 
the results are obtained from the simulations and the ac- 
tual uncertainty could be larger than these estimates. We 
exclude objects from the observed sample, if the loss is cal- 
culated to be larger than 90% in the profile fitting. In this 
paper, we use methodl as a fiducial line fitting method. 

3.2. Spec-z and Ha Line Flux 

The best-fit redshifts, which are referred to as spee- 
ds, are distributed from 1.2 < z < 1.6 with median of 
z = 1.408. In the upper left panel of Figure 5, the spec-z 
is compared to the photo-z. Although the spec-z gener- 
ally agrees with the photo-z with a 0.05, consistent with 
the uncertainty of the photo-z determination described in 
Section 2.2, there appears to be a systematic difference of 
~ 0.05 at 1.2 < z < 1.6. This may be caused by the sys- 
tematic uncertainty of the photometric zero-points. The 
presence of such systematic differences at some redshifts 
has been also reported by Ilbcrt et al. (2006). 

In Figure 5, we compare the fluxes of Ha and [Nil] 
emission lines obtained by the two different methods. 
Although the flux derived from the two methods gener- 
ally agree well for both emission lines, in some cases the 
flux measured by method2 is underestimated significantly, 
which is consistent with the results of the simulations, in- 
dicating the necessity of the correction for the OH mask 
effects. 

Comparisons of the observed Ha flux and the expected 
Ha flux are shown in Figure 6. The observed Ha fluxes 
are corrected for the fiber loss depending on their rso as 
we describe in Section 3.3. In the left panel, the expected 
Ha fiuxes calculated from the rest-frame UV luminosity, 
corrected for the nebular emission line using E{B — V), 
agree well with the observed Ha fluxes. In contrast, in the 
right panel, the expected fluxes by assuming the E{B~V) 
for the stellar component are systematically larger than 
the observed ones, implying that the dust extinction for 
the nebular emission is larger that for the stellar contin- 
uum. Figure 7 shows that the SFR derived from the ob- 
served Ha flux (extinction corrected for nebular compo- 
nent) agrees with that from the rest-frame UV luminosity 
density (extinction corrected for stellar component). 

Since the Ha emission line is superposed on the stellar 
absorption line, the correction for this effect is necessary 
in order to measure an accurate emission line fiux. Zahid, 
Kewley, & Bresolin (2011) fitted the stellar continuum of 
stacked spectra and found that the correction to H/? emis- 
sion is negligible in galaxies at z ^ 0.8. Since the continua 
of our sample are mostly too weak to be fitted, we esti- 
mate the absorption contribution from the SED models 
as follows. We calculate the equivalent width (EW) of 
Ha absorption line by assuming an exponentially decay- 
ing star formation history with r = 10 Myr and a stellar 
age of ~ 100 Myr; a typical case in our SED fitting for 
the observed sample. The resulting EW of — 4.2A is much 
smaller (- 2%) than the typical Ha EW of ^ 200A, and 
thus, we do not take the effect of stellar absorption into 
consideration. 



3.3. Aperture Effect 

A certain fraction of light outside the aperture of the 
FMOS fiber of 1"2 diameter is lost. Although the aperture 
effect on relative quantities such as a line ratio and the 
metallicity is small, the effect is critical to absolute quan- 
tities such as total flux and SFR. We estimate the cov- 
ering fraction for galaxies by using the NICMOS F160W 
archival images (Bouwens et al. 2010) convolved to the 
typical seeing sizcs(0."5, 0."75, 1"00, and 1"20 in FWHM) 
in the GOODS-N field. The object detection and l."2 
aperture photometry are carried out by using SExtractor. 
These objects are cross-matched with a K-selected cata- 
logue in GOODS-N/MODS (Kajisawa et al. 2010; Yuma 
et al. 2011). For objects brighter than Kg = 23.9 mag 
and more massive than M^, = 10^'^Mq at 1.2 < Zph < 1.6, 
we calculate the covering fraction, which is defined as 
FLUX_APER (in l."2 diameter) / FLUX.AUTO. The cov- 
ering fraction correlates very well with a half light radius 
deconvolved by the PSF (rso) in the F160W image; the 
scatter is very small (^ 4%). According to the simula- 
tions, the covering fractions are 0.20 — 0.80 depending on 
the seeing size; a typical value is ^ 0.5 when the seeing 
size is 0."75. 

In order to estimate fiber aperture losses for our sam- 
ple galaxies in the SXDS field, we firstly derive the ob- 
served based on the K-band WFCAM images by us- 
ing the SExtractor's FLUX_RADIUS. The intrinsic r^o of 
our sample is calculated from the observed rso from the 
K-band image, after deconvolution by a Gaussian PSF of 
0."80. The resulting rso distributes from 0."3 to 0"9 with 
the median value of 0."55 (4.6 kpc at z = 1.4). In order to 
examine the reliability of the estimation in K-band, these 
are compared to those obtained from higher-resolution 
images from WFC3/HST in the CANDELS survey held 
in a small part of the SXDS/UDS field. (Grogin ct al. 
2011; Koekemoer et al. 2011). The rso of our primary 
sample, described in Section 2.3, in the CANDELS FoV 
is measured from the F160W (the rest-frame wavelength 
is slightly different from that in K-band) images decon- 
volved by a PSF of 0."18. We found that they agree well 
with those derived from the K-band image within 0."05 
10% accuracy). 

The fiber loss for each object is determined from its 
rso. The seeing size is assumed to be 0."75 (FWHM) in 
i7-band, as expected for the seeing size of ^0"9 in R- 
band and the wavelength dependence cx A""'^ (e.g.. Glass 
1999). The median rso is 0."55 and the corresponding cov- 
ering fraction is 0.45. Line fluxes that we present in this 
paper are corrected for this aperture effect^. It should 
be noted that the positional error of the fibre allocation is 
<0."2 (Kimura et al. 2010) and the effect is small (<10%). 
Further, this error is random and does not cause a system- 
atic error. 

^ Strictly speaking, the aperture loss is already corrected by using 
the correction factor for a point source in the flux calibration. 
The correction factor should be different for a galaxy since the 
galaxy is generally extended. Although we need to modify the 
correction factor for a galaxy, the modification is less than 10% 
and we did not modify it. 
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Fig. 5. Comparison between spec-2 and photo-z {upper left). Comparisons of Ha flux {upper right), [N n] flux {lower left), and oxy- 
gen abundance {lower right) in two different spectral fitting methods. Open circles are objects whose S/N of [N ll] flux is less than 1.5. 



3.4. Possible AGNs 

Detections of Ha, [Nll]A6584, [O IIl]A5007, and H/3 lines 
allow us to make a line diagnosis to distinguish AGNs 
from star-forming (SF) galaxies (Baldwin et al. 1981). Ha, 
[Nil], [Olll], and H/3 lines are detected with S/N > 3 for 
7 objects. Ha and [Nil] lines are detected with S/N > 
3 but both [O III] and H/3 lines are not detected for 14 
objects. Both [OlIl] and H/3 lines are detected with S/N 
> 3 but Ha and [Nil] lines are not detected for 4 objects. 
The distribution of these objects on the BPT diagram is 
presented in Figure 8. The empirical criterion to separate 
the AGN and SF by Kauffmann et al. (2003) and the 
maximum theoretical line of starburst by Kewley et al. 
(2001) are plotted. Six objects are in the AGN region 
if the theoretical line is adopted. In addition to these 
objects, one object shows a line width larger than 1000 
km s~^, suggesting the presence of an AGN. Thus these 
seven objects are excluded from the sample hereafter. It is 
worth stressing that although several objects are located 
between the theoretical line and the empirical line, often 
referred to as a composite region, most are within the 
empirical line. 

In the right panel of Figure 8, we also show the stacked 



spectrum of all observed spectra excluding the seven AGN 
candidates. The stacking analysis shows that these galax- 
ies at z ^ 1.4 are located above the sequence of normal 
SF galaxies at z ^ 0.1. Similar trends have been reported 
at z - 2 (Erb et al. 2006; Hainline et al. 2009). By us- 
ing high S/N spectra of gravitationally lensed galaxies, 
Hainline et al. (2009) found that the galaxies at z 2 
show higher [O III] /[Oil] line ratios and thus higher ion- 
ization parameters compared to local SF galaxies. The 
systematic upward shift of our sample in the BPT diagram 
also suggests higher ionization parameters in galaxies at 
z ^ 1.4. For this reason, we use the maximum theoretical 
line by Kewley et al. (2001) as the selection of non-AGNs. 
However, we cannot completely rule out the possibility of 
contamination from low-luminosity AGNs. 

3.5. Metallicity Determination and Its Uncertainty 

Empirical indicators obtained from strong emission 
lines are used in order to derive the metallicity of high- 
redshift galaxies. The choice of the indicators depends on 
the targets' redshift. For galaxies at 1.2 < z < 1.6, H/3, 
[Qui], Ha, and [Nil] can be covered in the FMOS wave- 
length range. 



10 



K. Yabe et al. 




10"'° 10 
Observed Ha Flux (erg/s/cm^) 



10"'° 10" 
Observed Ha Flux (erg/s/cm^) 



Fig. 6. Comparisons of the observed Ha flux and the expected Ha flux. The expected Ho flux is calculated 
from the rest-frame UV luminosity corrected by E{B — V) for nebular {left) and stellar (right) components assum- 
ing the Calzetti et al. (2000) extinction law. Objects with >90% flux loss due to the OH-masks are indicated 
by red open circles. The AGN candidates in the BPT diagrams arc indicated by blue triangles (Section 3.4). 



A metallicity indicator with the line ratio of [N n] and 
Ha (N2 method) was calibrated by the following equation 
(Pettini & Pagel 2004): 



12 + log(0/H) = 8.90 + 0.57 x N2, 



(1) 



where N2 refers to log[/([Nn]A6584)//(Ha)]. The scat- 
ter around the relation is ~ 0.18 dex at 68% signifi- 
cance level (Pettini & Pagel 2004). It is also known 
that [Nn] tends to be saturated near and above solar 
metallicity. Pettini & Pagel (2004) found that 03N2 
(= log{([Oni]A5007/H/3)/([Nn]A6584/Ha)}) is more use- 
ful for higher metallicities. For 14 objects in our sam- 
ple, the [Nh], Ha, [Ohi] and H/3 lines are detected with 
S/N > 2. For these objects, we compared the metallic- 
ity derived by the N2 method to that from the 03N2 
method. We found that the metallicities derived from 
the 03N2 method are systematically lower by 0.07 dex 
than those from N2 method in the metallicity range of 
12-Hlog(0/H) = 8.4-8.8. The systematic difference 
of metallicity determination in these two methods by 
0.1 — 0.2 dex was reported at high-redshift (Erb et al. 
2006); this may be due to differing physical conditions 
in SF galaxies at z ^ 1.4 as compared to z ~ 0. 

In Figure 5, the metallicity derived from methodl is 
compared to that from methods. Although the metallic- 
ities derived from the two methods generally agree, the 
scatter is relatively large. 

3.6. Stacking Analysis 

More than half of the [N h] line fluxes are weak to mea- 
sure the [NH]/Ha line ratio and thus metallicity of in- 
dividual objects reliably. In order to cope with this low 
significance of the [N n] line and to reveal the 'average' na- 
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Fig. 7. Comparisons of the extinction corrected SFR ob- 
tained from the Ha luminosity and that from the 
rest-frame UV luminosity density. Objects with >90% 
flux loss due to the OH-masks are indicated by red 
open circles. The AGN candidates in the BPT di- 
agrams are indicated by blue triangles (Section 3.4). 
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Fig. 8. (left): BPT diagram. Upper and lower limits are 3(t values. The result from the stacking analysis is showrn by 
a black star. For the comparison, local SDSS galaxies are plotted as gray dots. The empirical criterion to separate the 
AGN and SF by Kauffmann et al. (2003) and the maximum theoretical line of starbursts by Kewley et al. (2001) are 
shown as dotted and dashed curves, respectively. (right): The stacked non-AGN spectrum for the H/3 and [Olll] region. 



ture of the sample galaxies, we separate the sample into 
three mass bins, stack the individual spectra, and derive 
the metallicity at the representative stellar mass in each 
bin. We perform the stacking analysis in two ways. One is 
stacking the best-fit spectra which are properly corrected 
for the OH-mask effect as described above, and the other 
is simply stacking the observed spectra. 

For the former method (stacking the model spectra), 
each model spectrum is made by a multi-Gaussian profile 
with the best-fit parameters from methodl (Section 3.1). 
These spectra are combined to give the average spectrum. 
The obtained spectra are again fitted by multi-Gaussian 
models to derive the line ratio and the metallicity. 

For the latter method (stacking the observed spectra), 
the observed spectra are averaged with weights based on 
the observed noise to minimize the effects of the OH-mask 
as follows. 



1 



(2) 



where Fi{\) and crj(A) are the de-redshifted observed flux 
density and noise, respectively. The noise is generally 
large at the OH masks. In the left panel of Figure 9, 
the spectra resulting from stacking the model spectra are 
shown. In the middle panel, the spectra are convolved 
to low resolution {R ~ 800 at the wavelength). In the 
right panel, the spectra resulting from stacking the ob- 
served spectra are shown. In the low resolution spectra, 
the spectra from the two different methods are roughly 
comparable. In Figure 9, it can be seen that the flux of 
[N II] increases with the increasing stellar mass. 

The uncertainties associated with the stacking methods 
are examined by stacking the sample in various ways. For 



both methods, we examine four different ways of combin- 
ing spectra: The weighted average, unweighted average, 
median, and flux normalized average. The typical uncer- 
tainties are 0.04 dcx (the largest mass bin) to 0.07 dex (the 
smallest mass bin). The uncertainties from the stacking 
method are almost comparable to or slightly larger than 
the observational errors. 

In the method of stacking observed spectra, some of 
the individual emission lines are affected by the OH- 
masks.This effect is simulated by stacking artificial spec- 
tra whose redshift is randomly distributed in the range of 
1.2 < z < 1.6. In this simulation, the line flux of Ha, the 
ratio of [Nil] /Ha, and the line width are assumed to be 
1.0 X 10~^^ erg s~^ cm~^, 0.20, and 200 km s~^, respec- 
tively, which are typical values for the detected objects. 
The number of the artificial spectra is ^ 20, which is com- 
parable to the number of the actual objects in each mass 
bin. We fit the realized spectrum and measure the line 
ratio and metallicity by method2 described in Section 3.1. 
This process is repeated 1000 times. The recovered line 
ratio agrees with the input value within ~ 20%, corre- 
sponding to ~ 0.1 dex in metallicity, with no systematic 
difference. 

In this paper, we use the former method, i.e., stack- 
ing the model spectra, as a fiducial stacking method. 
Hereafter, the stacked results are obtained by this method 
unless otherwise noted. 

3.7. The Mass-Metallicity Relation at z ^ 1.4 

The derived metallicity distribution as a function of 
stellar mass is presented in Figure 10. The stellar masses 
are re-calculated by the SED fitting in the same man- 
ner described in Section 2.2 but fixing the redshift to the 
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Fig. 9. Stacked spectra in three mass bins {red solid line) and the model fittings {blue dotted line). The profiles from 
the stacking analysis with best-fit models are in the left panels, those with best-fit models convolved to the LR spectra in 
the middle, and those with the "observed spectra" in the right panel. The mass ranges (the number of objects in each 
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from top to bottom. Vertical dotted linos indicate the positions of [Nll]A6548, Hq, and [Nll]A6584 from left to right. 




Fig. 10. Mctallicity as a function of stellar mass of our sample. The metallicity derived with spectral fitting with and without 
the effects of the OH masks are plotted in the left and right panel, respectively. The typical errors of stellar mass and metallicity 
are shown in the lower right corner of each panel. [Nil] lines with S/N > 3.0 and 1.5 < S/N < 3.0 arc indicated by filled and 
open symbols, respectively. [Nil] lines with S/N < 1.5 are plotted as upper limits with values corresponding to 1.5<T. Note 
that the error of the metallicity is derived from the flux error of Ho and [Nil] lines and does not include the uncertainty of the 
metallicity calibration. The results from stacking analysis are presented by filled stars. Thick solid lines are the second-order 
polynomials for the stacked results. Dashed lines and thin solid lines are the mass-metallicity relation at z ~ 0.1 (Tremonti et al. 
2004) and z ~ 2.2 (Erb et al. 2006). Horizontal dotted line indicates solar metallicity (12+log(0/H)=8.66; Asplund et al. 2005). 
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spectroscopic redshift obtained by the observations. All 
other physical quantities are based on the spec-z. In the 
left panel, metallicities are derived using methodl as we 
described in Section 3.6, while in the right panel, we use 
metallicities derived by meihod2. In both panels, the ob- 
jects for which the S/N of the [Nil] flux is less than 1.5 
arc plotted as upper limits with values corresponding to 
1.5(7. The mass-metallicity relations obtained at z ~ 0.1 
(Tremonti et al. 2004) and z - 2.2 (Erb et al. 2006) are 
plotted for comparison. Note that the stellar masses from 
these works are converted by a factor of 1.8 so that the 
IMF is consistent with that adopted in this paper. The 
metallicity at z ~ 2.2 is obtained by using the same N2 
indicator, and for that at z ~ 0.1 we use the recalculated 
values with the N2 indicator from Erb et al. (2006) for a 
fair comparison. In both panels of Figure 10, typical ob- 
servational errors of stellar mass and metallicity are shown 
at the lower-right corner. In these errors, the uncertainty 
of the metallicity calibration is not included. 

Since [N II] emission lines arc not detected significantly 
for ~ 40% of our sample, we measure the average metal- 
licity by dividing the sample by stellar mass and apply- 
ing stacking analysis including both detections and non- 
detections in the way presented in Section 3.6. We divide 
the sample into three mass bins, lO^ ^M© - 10^°-"Mo, 
1O1° "M0-1O1O-48M0, and IOIO ^^Mq - IO^ ooMq. The 
number of objects in each mass bin is ~20. The metal- 
licities derived from stacked spectra in three stellar mass 
bins are also presented as solid stars in the both pan- 
els of Figure 10. The errors are estimated by the boot- 
strap method with 1000 realizations. The resulting mass- 
metallicity relation obtained by our fiducial method (left 
panel) is systematically lower by 0.03, 0.04, and 0.05 dex, 
in the lowest, middle, and highest mass bin, respectively, 
than that from the observed spectra stacking method. 
The results are also summarized in Table 1. 

Figure 10 shows that most of individual objects and the 
results from the stacking analysis at z ^ 1.4 lie between 
the mass-metallicity relation at z ^ 0.1 and z 2.2. The 
stacked results are expressed by a second-order polyno- 
mial: 

12 + log(0/H) = -0.0724a;2 + 1.733a: - 1.67, (3) 

for the model spectra stacking method and 

12 + log(0/H) = -0.1082a;2 + 2.497a: - 5.71, (4) 

for the observed spectra stacking method, where x = 
log(M*/M0) in both equations. The resulting curves are 
also presented in both panels of Figure 10. 
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Fig. 11. Intrinsic scatter of the mass-metallicity relation 
at 2 ~ 1.4 {filled stars). Strictly speaking, this is a 
lower limit of the scatter. Error bars are based on 
the bootstrap method. Scatters at 2 ~ 0.1 by Tremonti 
et al. (2004) are indicated by filled circles and theo- 
retical predictions at 2 = 0, 1, 2, and 3 by Dave, 
Finlator, & Oppenheimer (2011) are indicated by solid lines. 

In section 3.3, we found that the typical covering frac- 
tion of FMOS fiber is ^ 0.5. The covering fraction by 
the fibre aperture affects the estimation of metallicity if 
there is a steep radial metallicity gradient; a small cover- 
ing fraction leads to an overestimation of the metallicity 
if the metallicity gradient is negative. Tremonti et al. 
(2004) showed that the aperture effects on the metallicity 
is ^ 0.1 dex by using the SDSS galaxies with the typical 
covering fraction of 0.24. In the local universe, the dif- 
ference between the metallicity measured in the aperture 
and the integrated metallicity is ^0.1 dex for the covering 
fraction of > 0.2 (Kewley et al. 2005). Since the typical 
covering fraction of our sample is 0.5, the metallicity 
gradient is unlikely to affect the derived metallicity so 
much unless the metallicity gradient drastically changes 
at high- redshift. Recently, steeper metallicity gradients 
than local galaxies have been reported at z ^ 2 (Jones 
et al. 2010; Yuan et al. 2011), while the flatter or even 
"positive" gradients have also been reported (Cresci et al. 
2010). The metallicity gradient at high redshift remains 
unclear due to the small size of currently available sam- 
ples. 

The intrinsic scatter of the mass-metallicity relation has 
been reported at z 0.1 (Tremonti et al. 2004; Yates 
et al. 2011) and also reproduced by cosmological simu- 
lations (Dave, Finlator, & Oppenheimer 2011). We try 
to constrain the scatter of the mass-metallicity relation at 
z ^ 1.4, though the sample size may not be large enough 
yet. Figure 10 shows that the metallicity distribution has 
a larger scatter as compared with the typical observational 
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Fig. 12. Dependence of SFR on the mass-metallicity relation. In the left panel, the sample is divided into three groups 
with their SFRs derived from Ha luminosity corrected for dust extinction: Circles for detection and arrows for non-detec- 
tion for [Nil] with SFR < 53.0 Moyr"! {small red), 53.0 < SFR < 85.0 A/oyr^i {middle green), SFR > 85.0 Moyr"! 
{large blue). In the right panel, the sample is divided by their SFRs derived from the rest-frame UV luminosity density cor- 
rected for dust extinction, with SFR < 40.0 Mqyt-^ {small red), 40.0 < SFR < 80.0 MQyr'^ {middle green), SFR > 80.0 
M0yr~^ {large blue). The filled stars show the stacked results in the higher SFR group, while the open stars show those in 
the lower SFR group in each mass bin. Solid curves show the mass-metallicity relations described by Eq. 3 in Section 3.7. 



error. So the scatter is defined as the RMS of the devia- 
tion from the best-fit mass-metallicity relation (Eq. 3) in 
each mass bin, after subtracting the contribution from the 
averaged observational error of individual objects. The re- 
sulting standard deviations are 0.138, 0.113, and 0.057 
dex in mass bins of log(Af*/Afo) = 9.98, 10.28, 10.65, re- 
spectively. Since we take the metallicities from upper lim- 
its for many objects, the intrinsic scatter should be larger 
than these values. Figure 11 shows the intrinsic scatter 
against stellar mass. The scatters of the mass-metallicity 
relation at z ~ 1.4 are comparable to those obtained at 
z ~ 0.1 by Tremonti et al. (2004); they decrease with in- 
creasing stellar mass of galaxies if we see the face values, 
though the small scatter at the massive part may be partly 
due to the saturation effect of the N2 indicator. 

3.8. Dependence of SFR and Size on the Mass- 
Metallicity Relation 

We examine the dependency of SFR on the mass- 
metallicity relation for our sample at 2; ^ 1.4. In Figure 
12, we present the metallicity as a function of stellar mass 
for individual objects dividing into three groups accord- 
ing to their intrinsic SFRs. In the left panel, the SFR 
is derived from the observed Ha luminosity corrected for 
extinction by using E{B — described in Section 2.2, 

while in the right panel, the SFR is derived from the rest- 
frame UV luminosity density also corrected for extinction. 
They show no clear trend that galaxies with higher SFRs 
show lower metallicities, rather, they show a correlation 
between the SFR and the stellar mass, i.e., galaxies with 



larger SFRs have larger stellar masses. In order to test 
this further, we split the sample in each mass bin into two 
groups by the median SFR. In each mass bin and SFR 
group, the best-fit model spectra obtained by methodl 
are combined. The results are presented in Figure 12 as 
filled and open stars for higher and lower SFR groups, re- 
spectively, as well as in Table 2. For both SFRs derived 
from Ha and UV, the resulting metallicities in higher SFR 
groups are lower than those in lower SFR groups in all 
mass bins, though the difference at the most massive bin 
in the right panel is very small. At A/, = lO^^Af© in the 
Figure 12, the gradient of metallicity against SFR from 
Ha is A(12-hlog(0/H))/A(log(SFR)) - -0.6 dex, which 
is smaller (steeper) by a factor of ~ 2 than that of ^ —0.3 
dex at 2; ~ 0.1 (Mannucci et al. 2010). The metallicity 
gradient against SFR from UV is ^ —0.4 dex. which is 
almost similar to that at z 0.1. 

Figure 13 shows that the dependence of specific SFR, 
derived from Ha (left panel) and the rest-frame UV (right 
panel), on the mass-metallicity relation is similar to that 
of the SFR (Figure 12), but the trend is relatively unclear. 
This is probably because the correlation between the stel- 
lar mass and the SFR cancels out the dependence on the 
metallicity. 

We also examine the dependency of galaxy size (half 
light radius; r^o) described in Section 3.3 on the mass- 
metallicity relation. In Figure 14 and Table 3, in the 
same manner as the SFR, we present the metallicity as 
a function of stellar mass dividing the sample by the size. 
Both individual and the results from stacking show that 
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Fig. 13. Dependence of specific SFR (sSFR) on the mass-metallicity relation. In the left panel, the sample is divided into three 
groups according to sSFRs derived from Ho luminosity corrected for dust extinction: Circles for detection and arrows for non-detec- 
tion for [Nil] with sSFR < 2.5 Gyr'^ (small red), 2.5 Gyr"! < sSFR < 4.0 Gyr"! {middle green), sSFR > 4.0 Gyr"! (large 
blue). In the right panel, the sample is divided according to sSFRs derived from the rest-frame UV luminosity density cor- 
rected for dust extinction, with sSFR < 2.3 Gyr^^ (small red), 2.3 Gyr'^ < sSFR < 3.6 Gyr^^ (middle green), sSFR > 3.6 
Gyr~^ (large blue). The filled stars show the stacked results in the larger sSFR groups, while the open stars show those in 
the smaller sSFR groups in each mass bin. Solid curves show the mass-metallicity relations described by Eq. 3 in Section 3.7. 



Table 2. Metallicity from the stacked spectra in three stellar mass bins and two SFR groups 



log(A/*/M0) 


gpj^mcd (jj^) 


Spj^med (UV) 




12-Mog(0/H) 






Low SFR (Ha) 


High SFR (Ha) Low SFR (UV) 


High SFR (UV) 




56.7 


25.8 


8.493±0.049 


8.286±0.071 8.420±0.067 


8.337±0.087 


10.28i^-20 


55.9 


54.0 


8.511±0.050 


8.468±0.056 8.571±0.057 


8.420±0.040 


10.65t°?| 


112.9 


91.5 


8.614±0.020 


8.556±0.042 8.578±0.034 


8.572±0.020 



galaxies with larger sizes tend to have lower metallicities. 
In the stacking analysis, the difference is larger as the 
stellar mass decreases. At AI^, = IO^'^Mq, the gradient of 

metallicity with rjo is A(12 Iog(0/H))/A(r5o) 0.1 

dex kpc~^, which is comparable to ~ —0.1 dex kpc~^ ob- 
tained at 2; ~ 0.1 by Ellison et al. (2008a). Since our sur- 
vey covers the CANDELS field where deep HST images 
arc available, we can also consider the morphological in- 
formation for a part of our sample. The number of such 
galaxies, however, is too small to be considered statisti- 
cally. The relation between the metallicity and morphol- 
ogy will be discussed in the future works. 

We further examine the dependency of other physical 
quantities, such as colors, surface mass density, color ex- 
cess. No clear trend, however, can be seen. It is not 
certain whether there are truly no trends for such param- 
eters because the sample size is still not large enough. A 
larger sample from the further survey with FMOS in the 
future is desirable in order to study these questions. 



3.9. Effects of the Sample Selection on the Stellar-Mass 
Relation 

As mentioned in Section 2.3, our targeted sample is 
somewhat biased toward the larger SFR end of the galaxy 
sequence in the Af* — SFR diagram as compared to the 
primary sample. According to the dependency of SFR on 
the mass-metallicity relation at z ~ 1.4 described above, 
this selection may lead to a bias toward the smaller metal- 
licity. The expected metallicity of the primary sample 
at a given stellar mass is examined by assuming the de- 
pendency of the SFR on the mass-metallicity relation de- 
scribed in the previous section. The metallicity is ex- 
pressed as a linear function of intrinsic SFR from UV in 
each mass bin and expected metallicities of primary sam- 
ple are estimated from their SFRs. The average metallic- 
ity weighted with the Ha flux is still close to the metal- 
licity derived in this study within ^ 0.02 dex. 

Our sample is also slightly biased toward a larger size as 
compared with the primary sample. We also examine the 
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Fig. 14. Dependence of half light radius (?'5o) on the 
mass-metallicity relation. The sample is divided into three 
groups by their half light radii with rso < 4.38 kpc {small 
red), 4.38 kpc < r$Q < 5.30 kpc {middle green), rso > 5.30kpc 
{large blue). The filled stars show the stacked results in the 
larger sSFR (rso) groups, while the open stars show those in 
the smaller rso groups in each mass bin. Solid curves show the 
mass-metallicity relations described by Eq. 3 in Section 3.7. 



Table 3. Metallicity from the 
in three stellar mass bins and 



stacked spectra 
two rso groups 



log{AU/MQ) 




12+log(0/H) 






Small Large rso 




4.4 


8.477±0.041 8.365±0.053 


10.28;°?° 


4.7 


8.513±0.050 8.479±0.057 


10.65;EJ;?? 


5.6 


8.562±0.017 8.577±0.035 



effects on the metallicity by assuming the dependency of 
size on the mass-metallicity relation in the same manner 
described above. The resulting metallicity is also close 
to the observed metallicities within ~ 0.03 dex. These 
effects of both SFR and size on the metallicity are smaller 
than the other uncertainties. Hence, we conclude that 
the effects of the selection bias on the mass-metallicity 
relation are relatively small. 

The color excess of our sample distributed widely from 
0.0 mag to ~ 0.6 mag with median value of ^0.23 mag. 
The sample is slightly biased toward a smaller color excess 
as compared with the primary sample (from 0.0 mag to 
'^1.0 mag with the median value is 0.26 mag). However, 
since there is no clear dependence of color excess on the 
mass-metallicity relation, the selection effect on the mass- 
metallicity relation is small. 



4-.1. Consecutive Evolution of the Mass-Metallicity 
Relation from z ~ 3 to z ^ 

In Figure 15, the mass-metallicity relation obtained in 
this study is shown with those at z ~ 0.1 (Tremonti et al. 
2004), at z - 0.8 (Zahid, Kewley, & Bresolin 2011), at 
z 2.2 (Erb et al. 2006), and at z ~ 3.1 (Mannucci et al. 
2009). Note again that the stellar mass and the metallicity 
of these samples are converted so that the IMF and metal- 
licity calibration are consistent with those we adopted. 
For the conversion of the metallicity, we use prescriptions 
by Kewley & Elhson (2008) for the results at z - 0.8 and 
by Nagao, Maiolino, & Marconi (2006) for z 3.1. 

In Figure 15, it can be seen that our result at z ^ 1.4 
is located between that at z 0.8 and z ^ 2.2 and the 
mass-metallicity relation evolves consecutively from z 3 
to z 0.1. There appears to be no clear downsizing evo- 
lution from z ~ 3 to z 0.8. It is interesting to note 
that there seems to be a downsizing-like evolution from 
z = 0.8 — 1.4 to z ~ 0.1. The metallicity difference be- 
tween z ^ 1.4 and z ~ 0.1 is ~ 0.2 dex at ~ W^^Mq, 
^ 0.1 dex at 1O^°'^M0, and almost comparable at 

M, ~ lO^^Af0 if we extrapolate our mass-metallicity rela- 
tion (Eq. 3). However, here we recall the saturation effect 
of the N2 indicator. Actually, the recalculated result of 
SDSS galaxies from Erb et al. (2006) shows a flattened 
MZ-relation at ^ 1 x IO^^Mq for near-solar metallic- 
ity. Since the metallicities of our sample at z 1.4 are 
mostly sub-solar, they are not expected to suffer from se- 
vere saturation effects. The flat feature at the high mass 
end at z 0.1, however, may be due to the saturation 
effect. 

Since the mass metallicity relation evolves without 
changing its shape so much except for the massive part at 
z ~ 0.1, we trace the chemical evolution at M* = 1O^°M0. 
In the upper panel of Figure 16, the metallicity at = 
1O^°M0 from the mass-metallicity relation at each red- 
shift is plotted as a function of redshift. It clearly shows 
that the mean metallicity increases with decreasing red- 
shift. The cosmological evolution is well reproduced with 
12 + log(0 /H) = 8.69 - 0.086(1 + z)!-^. In the lower panel, 
the metallicity is also plotted against the cosmic age. The 
metallicity evolution rate changes rapidly at z ^ 2 and 
seems to be saturated after that. Interestingly, the end 
of the rapid growth of metallicity roughly corresponds to 
the peak epoch of the cosmic star-formation history. It 
is worth noting here that the sample selections for these 
mass-metallicity relations are not necessarily the same. 
Although the mass-metallicity relation may be affected 
by the sample selection, pursuing the effect at high red- 
shift is still challenging and will be studied in the future 
works. 

The mass-metallicity relation of our sample is also com- 
pared to that obtained from recent cosmological hydro- 
dynamic simulations by Dave, Finlator, & Oppenhcimer 
(2011). They implemented various outflow process and 
found that the momentum-conserving winds, in which 
both the outflow rate and the wind velocity depend on 
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Fig. 15. Comparisons to the mass-metallicity relations at 
z~0.1 {thin solid) by Tremonti et al. (2004), 2~0.8 (dashed) 
by Zahid, Kewley, & Bresolin (2011), 2 ~ 2.2 {dotted- dashed) 
by Erb et al. (2006), and z ~ 3.1 {dotted) by Mannucci 
et al. (2009). The ranges of the relations show the ob- 
served range of the stellar masses. The stellar masses and 
the metallicitics of other samples are converted so that the 
IMF and metallicity calibration are consistent with those we 
adopted. Horizontal dotted line indicates solar metallicity. 



the velocity dispersion of the galaxy, i.e., galaxy mass, 
could reproduce the overall shapes of the mass-metallicity 
relations at z~0.1, z~2, and z 3 well (also shown in 
Figure 16). In Figure 17, we show the mass-metallicity 
relation obtained at 2; ^ 1.4 as well as those in the simu- 
lations with four wind models: No winds {nw), constant 
winds (cw), slow winds (sw), and momentum-conserving 
winds (vzw). We take the average of results at z ~ 1 and 
z ~ 2 for comparison at z ~ 1.4. Figure 17 shows that our 
mass-metallicity relation agrees well with the vzw model, 
although it also agrees with the cw model within the er- 
ror bars. This agreement between our mass-metallicity 
relation and the theoretical vzw or cw models may possi- 
bly indicate the existence of strong and mass-dependent 
winds in galaxies at z = 1 — 2. 

In Figure 11, the scatter of the mass-metallicity relation 
at z 1.4 is compared with the predictions at z = — 3 with 
the vzw model by Dave, Finlator, & Oppenheimer (2011). 
While the model shows the evolution of the scatter from 
z = l — 2 to z = 0, our result shows no clear evolution from 
z ~ 1.4 to z 0.1. The estimated scatter in this study, 
however, is a lower limit. If the scatter aX z ^ 1.4 is sig- 
nificantly larger than that shown in Figure 11, the scatter 
was larger at the redshift than that seen at z ~ 0.1; the 
scatter of the mass-metallicity relation may decrease as 
the redshift decreases. This trend is opposite to the simu- 
lation by Dave, Finlator, & Oppenheimer (2011). Further 
deep observations with a larger sample are desirable. 




4 6 8 10 

Cosmic Age (Gyr) 

Fig. 16. Evolution of the mean metallicity at M* = IO^^Mq 
against redshift {upper) and cosmic age {lower). The metal- 
licitics at z ~ 0.1 (Tremonti et al. 2004), z 0.8 (Zahid, 
Kewley, & Bresolin 2011), 2 ~ 1.4 (this work), 2 ~ 2.2 (Erb 
et al. 2006), and 2 ~ 3.1 (Mannucci et al. 2009) are used. 
In the upper panel, the linear regression line is also pre- 
sented as a solid line: 12 -|- log(0/H) = 8.69 - 0.086(1 -|- 2)1-3. 
In both panels, the vzw models by Dave, Finlator, & 
Oppenheimer (2011) are also shown as dashed lines. 



4-2. A Possible Scenario for the Dependence of Various 
Parameters on the Mass-Metallicity Relation 

In Section 3.8, we presented the dependency of the 
(specific) SFR and galaxy size on the mass-metallicity 
relation. It is suggested that the SFR or specific SFR 
is a second parameter which causes the scatter of the 
mass-metallicity relation both observationally (Ellison 
et al. 2008a; Mannucci et al. 2010; Lara-Lopez et al. 
2010; Yates et al. 2011) and theoretically (Dave, Finlator, 
& Oppenheimer 2011; Yates et al. 2011) at z 0. By 
using SDSS galaxies at z '-.^ 0.1, Mannucci et al. (2010) 
found that galaxies with larger SFR tend to show lower 
metallicity, and the scatters are significantly reduced by 
their proposed relation introducing SFR as the second pa- 
rameter, i.e., the fundamental metallicity relation (FMR). 
They also found that the mass-metallicity relations ob- 
tained at z = 0.5 — 2.5 are consistent with the FMR at 
z 0.1; the evolution of the mass-metallicity relation from 
z ~ 2.5 to z ^ may be explained by the difference of SFR 
of the samples. In Figure 18, we plot our data points on 
the fundamental metallicity relation by Mannucci et al. 
(2010), i.e., metallicity against log(M*)-^log(SFR), where 
/i is a projection parameter. Mannucci et al. (2010) found 
that ^ = 0.32 gives the minimum intrinsic dispersion of 
metallicity. The metallicity scatter of our sample is not 
reduced radically by the projection. The average metallic- 
itics in all mass bins are close to the FMR, but are shifted 
slightly; ji = 0.22 gives the best fit. It should be noted. 
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Fig. 17. Comparison to the theoretical models by Dave, 
Finlator, & Oppenheimer (2011). The theoretical pre- 
dictions for 2 ~ 1.4 are obtained by averaging those 
for z = 1.0 and z = 2.0. No winds (nw), Constant 
winds (cw), Slow winds (sw), and Momentum-conserving 
winds (vzw) are indicated by gray, green, orange, and 
cyan, respectively. The observed mass-metallicity rela- 
tion is the same as that in the left panel of Figure 10. 



however, that SFRs in our sample are mostly larger than 
those in the FMR at z 0.1. 

Various interpretations for the large scatter of the mass- 
metallicity relation and the dependence of second param- 
eters have been discussed. One of the possible interpreta- 
tion is the dilution effect of inf ailing gas. The infall of the 
gas from IGM through cold (smooth) accretion or galaxy 
mergers is considered to be important in galaxy growth 
(e.g., Bouche et al. 2010; Hopkins et al. 2006). When 
pristine gas falls into a galaxy, preexisting enriched gas is 
diluted by the fresh gas, and at the same time, the new 
star-formation is likely to occur, raising the SFR. 

The size dependence can also be explained by the infall 
of gas. It is suggested that a galaxy merger event not only 
can supply fresh gas but also may be able to increase the 
effective radius of the galaxy (Bezanson et al. 2009; van 
Dokkum et al. 2010). A correlation between the signal 
of merger events and the metallicity has been suggested 
(e.g., Elhson et al. 2008b; Rupke et al. 2010; Sol Alonso 
et al. 2010). According to their results, interacting or 
merging galaxies tend to show lower metallicity than iso- 
lated galaxies. In our sample, we found that 20 galaxies 
are closely paired according to a pairing analysis (e.g., 
Patton et al. 2002) with requirements that (1) the pro- 
jected separation is < 30 kpc (2) the redshift separation 
is < 0.15, which corresponds to 3(t of photo-z uncertainty. 
For these objects, however, no clear trends that the inter- 
acting galaxies tend to be metal-poor can be found. 




9.5 10 10.5 
log(M,)-0.321og(SFR) 

Fig. 18. The metallicity distribution of our sample against 
log(M») — 0.321og(SFR), i.e., the fundamental metallicity re- 
lation {solid line) proposed by Mannucci et al. (2010). Data 
points are based on those in the left panel of Figure 10. 



5. Conclusions and Summary 

We present the first results obtained from near-infrared 
spectroscopic observations of star-forming galaxies at z ^ 
1.4 with FMOS on the Subaru Telescope. We observed 
K-band selected galaxies in the SXDS/UDS fields with 
K < 23.9 mag, 1.2 < Zph < 1.6, M, > IO^-^Mq, and ex- 



pected F(Hq!)> 10~ erg s~ cm^ . Ha emission lines 
of 71 objects are detected significantly. For these ob- 
jects, excluding possible AGNs identified from the BPT 
diagram, gas-phase metallicities and upper limits are ob- 
tained from [Nil] /Ha line ratio by carefully considering 
the effects of the OH-masks. We separate the sample into 
three stellar mass bins and stack the spectra. We obtain 
a mass-metallicity relation at z ~ 1.4 from the stacking 
analysis. The mass-metallicity relation is located between 
those at z ~ 0.8 and z ^ 2.2. We tried to constrain an in- 
trinsic scatter of the mass-metallicity relation and found 
that the scatter is ^0.1 dcx and is comparable to that 
at z '--^ 0.1 (Trcmonti et al. 2004). The scatter increases 
as stellar mass decreases if we take the intrinsic scatters 
at face-value. The scatters of the mass-metallicity rela- 
tion at z ^ 1.4, however, should be lower limits; this im- 
plies that the scatter may be larger at higher redshifts. 
We found that the deviation from the mass-metallicity 
relation depends on the SFR and the half light radius: 
Galaxies with higher SFR and larger half light radii show 
lower metallicities at a given stellar mass. One possible 
scenario for the trends is the infall of pristine gas accreted 
from IGM or through merger events. Our results do not 
show clear FMR as proposed by Mannucci et al. (2010). 
Trends of the dependence of the SFR and the size on the 
mass-metallicity relation are, however, still not so clear. A 
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larger sample from the further survey with FMOS in the 
future may be able to reveal clearer trends, not only the 
dependence of the SFR and the size but also that of other 
parameters. The compilation of the mass-metallicity rela- 
tions at z 3 to z ~ 0.1 shows that they evolve smoothly 
from z 3 to 2 without changing its shape so much 
(except for the massive part at z ^ 0.1). The metallic- 
ity at = 10^°Mq on the mass-metallicity relation in- 
creases with decreasing redshift and can be described as 
12 -f- log(0/H) = 8.69 - 0.086(1 -I- z)^-^. The metallicity 
evolution rate was the highest at the cosmic age of ^ 3 
Gyr which was before redshift of 2. However, the result 
may be influenced by sample selection and further studies 
with large samples at high redshifts are desirable. 
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